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Abstract 

We compute complete tree level matrix elements for gg, qq — > hhW^W~ . We 
analyze the irreducible backgrounds to top signal at the Tevatron and at the 
LHC Their contribution to the total cross section is about 5% at the LHC, 
due to single resonant channels. Several distributions with contributions 
from signal and backgrounds are presented. 
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Introduction 



With the discovery of the top quark at the Tevatron [|l| by the CDF and DO col- 
laborations, all fermions belonging to the three generations which, according to the 
LEP measurement of the invisible width of the Z, possess a light neutrino have been 
observed. 

Using 67 pb~^ of data CDF measured the top mass to be Mf = 176 ± 8(stat) ± 
lO(syst) GeV and the production cross section as au = Q-Slil'^ ph. The larger sample of 
top events which has been subsequently collected will soon allow a better determination 
of the parameters and properties of the newly discovered particle. Already at this year 
winter conferences new preliminary results have been presented, based on about 110 
pb~^ of data 0. The most precise measurement of the top mass, extracted from the 
lepton plus three or more jets sample, is Mt = 175.6 ±5.7(stat) ± 7.1(syst) GeV. The 
corresponding result for the production cross section is au = 7.51} g pb. 

The mass of the top is a fundamental parameter of the Standard Model and plays 
a crucial role in radiative corrections to electroweak observables. In fact bounds on the 
top mass can be obtained from a comparison of high precision measurements at LEP 
and SLC, together with information from neutrino scattering, with the best theoretical 
predictions. Within the Standard Model the top mass must lie in the range Mt = 
178 ±8^20 GeV 0. The agreement between this indirect determination of Mt with the 
value which is directly measured at Fermilab is a remarkable success of the model. Any 
improvement in the precision of the top mass measurement, particularly when the mass 
of the W boson will be measured to a precision of about 50 MeV at LEP 2, will further 
test the electroweak theory, significantly reducing the allowed range for the mass of the 
Higgs boson. 

It is obvious that the discovery of the top is only the first step. All details of the 
production process and decay will have to be thoroughly examined. The measurement 
of the production cross section and of the distributions of different kinematic variables, 
like the pt of the top, the angular distribution of the decay products or the character- 
istics of additional gluon jets, will challenge our understanding of perturbative QCD. 
The determination of the top decay channels will test whether the new heavy quark be- 
haves as predicted by the Standard Model and might provide a window on new physics. 
These studies will be continued at Fermilab after the construction of the Main Injector 
and will be further refined at the Large Hadron Collider (LHC), where the cross section 
is much larger than at Tevatron energy and at the Next Linear Collider (NLC), in the 
clean environment of an e^e~ machine. 

In order to achieve a complete understanding of top production, a theoretical effort 
matching the advances on the experimental side is required. The production cross 
section has long been computed at next to Leading Order (NLO) in QCD [Q. More 
recently the contribution of soft gluons have been resummed to all order |^ using 
different techniques. In ref. the resummation procedure has been applied to the 
inclusive transverse momentum and rapidity distributions of the top. Several kinematic 
distributions have been computed at NLO in ref. [|^. All these studies have treated 
the top quark as a stable particle, separating the production of a it pair from the 
independent decay of the two heavy quarks, in the spirit of the so-called narrow width 
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approximation (NWA). 

The overall uncertainty in the calculation of the total cross section due to different 
choices of parton distribution functions and of renormalization and factorization scales 
is of the order of 20%. The difference between pure NLO calculations and calculations 
which also resum soft gluon effects is also of order 20% at the Tevatron and as large 
as a factor of two at the LHC These uncertainties however affect in a similar way all 
the contributions to the total cross section, including irreducible backgrounds whose 
relevance remains to be assessed. 

In this letter we examine the production of the hhW^W~ final state at the Tevatron 
and at the LHC, taking into account the full set of tree level diagrams. Our approach 
allows us to study the effects of the finite width of the top and of the irreducible 
background to it production, including their interference. Correlations between the 
decays of the two t-quarks are automatically included. 

A few representative diagrams are shown in figs. 1 and 2. In addition to the standard 
mechanism which is 0{a'^a1) we have studied the contribution of the 0{a'^) diagrams 
in which an initial pair of light quarks annihilates to a photon or a Z boson. Many 
of the diagrams we are studying (3,4,5 in fig. 1; 2 in fig. 2), while contributing to the 
irreducible background to ti production, effectively represent the production of a single 
top and must be taken into account in any attempt to study this rarer top production 
mechanism which is of interest in his own right. The simpler process qq — > W* 
tb, for instance, has been discussed in ref as a promising candidate for a precise 
measurement of the Vtb element of the CKM mixing matrix. In ref. |^ the subset of 
diagrams (1,2 plus crossed in fig. 1 and 1 in fig. 2) which describe the production of a 
it pair qq, gg — >• it followed by the two decays t — ^ W^b, i — > W~b has been studied 
using helicity amplitude methods. Strictly speaking, since the three diagrams discussed 
in ref. are not by themselves gauge invariant, the results obtained in this fashion 
are doubtful. It can however be argued that using an appropriate gauge, and under 
experimental conditions optimized for the observation of on-shell top pairs, the error 
will be of order Tt/Mf. Here we extend the helicity amplitude approach to the full 
gauge invariant set of diagrams which are required to describe bbW^W~ production. 

Within the present range of experimentally allowed top mass values the top 
lifetime is of the order of 10~^^ sec and the top decays to bW^ before hadronizing. 
As a consequence, there might be significant finite width and irreducible background 
effects at the percent level. Indeed we have previously found corrections of this order of 
magnitude in e+e~ bbW^W~ at the NLC [1^. It is therefore important to evaluate 
these effects and, when relevant, fully include them in the theoretical analysis. The 
determination of the top mass and of the production cross section are clearly influenced 
by irreducible backgrounds and one cannot rely on the NWA without checking the 
size of the corrections to this approximation. Furthermore, signal and background 
usually have different distributions in phase space and the spectra of some interesting 
kinematical observables may be more sensitive to these corrections than the total cross 
section. Alternatively, these differences may be useful for a more effective separation 
of the top signal from non resonant backgrounds. 

We have made no approximation concerning the mass of the 6, taking it into full 
account both in phase space and into the computation of the matrix elements. The b 
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mass also regulates the collinear singularities that would be present in some of the non 
resonant diagrams, for instance diagram 3, 4, 6 and 7 in fig. 1, if we had treated the b 
quark as massless. 



Calculation 

The matrix elements for gg W~^W~bb and for qq — > W~^W~bb have been computed 
using helicity amplitude methods. For the first process the method of ref. has 



been used. The amplitude for the second reaction has been computed with the method 
described in ref. . The expression of the corresponding Feynman diagrams has been 



generated with the help of PHACT (Program for Helicity Amplitude Calculations with 
Tau matrices) |T^. We have checked our results for gauge and BRST invariance [|14 . 



The full set of diagrams has been divided into four subsets, which were named 
after the corresponding structure of resonant enhancements as double resonant, t single 
resonant, t single resonant, and non resonant. 

Some representative diagrams for gg — > W^W~bb are shown in fig. 1. The full set 
includes 39 0{a'^a1) diagrams including 8 diagrams with a virtual Higgs. These latter 
are suppressed by the small coupling of the Higgs boson to the b and can be neglected. 
Diagrams 1 and 2 are double resonant. They correspond in narrow width approximation 
to the production and decay of a ti pair. Integrating over the whole phase space the 
amplitude squared obtained from just these two diagrams, one reproduces in NWA the 
cross section for gg it. The branching ratio of t (i) to W'^b {W~b) is in fact one, as 
we take off-diagonal CKM matrix elements to be zero. 

Diagrams 3, 4 and 5 are examples of diagrams which are single resonant in the t 
channel. They correspond in narrow width approximation to the production of a t- 
quark together with a bW~ pair followed by the decay t bW'^. All diagrams in 
the same class can be obtained permuting gluon lines and connecting the negatively 
charged W in all possible positions, while leaving the next to the outgoing fermion. 
The single resonant diagrams corresponding to the production of a t particle can be 
easily obtained in a similar way. The remaining diagrams in fig. 1, numbered 6, 7 and 
8, are examples of non resonant diagrams. 

Selected diagrams describing qq W^W~bb are shown in fig. 2. The full set 
includes 16 Oi^a^al) diagrams including 2 diagrams with a virtual Higgs and 62 0{a^) 
diagrams, of which 11 include a Higgs. In the former case the diagrams with a Higgs 
propagator are suppressed by the Hbb couphng. In the latter case the Higgs can be 
connected with a large coupling to the intermediate Z. However the contribution of 
these diagrams can be relevant only if the Higgs mass is larger than 2Mw and the 
Higgs propagator can go on mass-shell. To simplify our argument, we have neglected 
all diagrams which include a Higgs particle, limiting our discussion to the diagrams 
which survive in the limit of an infinitely heavy Higgs. In case a Higgs is discovered, 
the complete set of diagrams would have to be considered. For a similar discussion at 
the NLC see ||10| , p!5[| . In analogy with the gg W^W~bb case, the full set of diagrams 



can be subdivided in a double resonant (1 of fig. 2) subset, two single resonant subsets 
corresponding to t production (2 of fig. 2) and i production and a non resonant subset 
(3,4,5,6). 
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We have found it convenient to separate the different contributions to the cross 
section and to integrate them separately, adapting the choice of phase space variables 
to the peaking structure which correspond to the different sets of intermediate particles 
which can go on mass shell. We have identified a double resonant contribution gen- 
erated by the modulus square of the corresponding diagrams and two single resonant 
contributions which include the modulus square of the single resonant diagrams plus 
the interference between single and double resonant diagrams. All remaining terms, 
the modulus square of non resonant diagrams plus the interference between single res- 
onant and non-resonant diagrams, give the last, non-resonant, contribution. For both 
reactions the four terms have been numerically integrated using Vegas [1^ and the 
resulting distributions have been summed in the end. 

The relative size of the different contributions is obviously dependent on the gauge 
and only the sum of all terms is gauge independent. In a given gauge however one can 
meaningfully discuss the numerical significance of the various expressions. In particular, 
in the physical, purely transverse gauge we have employed, there is a clear hierarchy 
between the four contributions in which we have divided the full cross section. The 
double resonant contribution is always much larger than the single resonant ones, which 
in turn dominate the non resonant term. 

We have verified numerically that in NWA the total cross sections calculated from 
the double resonant subsets reproduce the results obtained for gg it and qq it for 
on shell top-quarks. 

We have used the MRS(D_) Jl^ set of parton distribution functions throughout this 



paper. For we have used the one-loop expression with Q = M^. A different choice 
of parton distribution functions or of the scale at which is evaluated would slightly 
alter our results for the various cross section. However we are here mainly interested in 
the relative weight of the different contributions and their ratio is essentially insensitive 
to these choices. 



Results 

Our results for the total cross section are presented in table 1 and table 2. For the 
Tevatron we have studied only the energy available at present, ^/s = 1.8 TeV, neglecting 
possible improvements that have been recently discussed in the literature. The LHC 
will probably begin its operation with a center of mass energy of 10 TeV which will 
be raised after a few years of running to y/s = 14 TeV. Both energies have thus been 
examined. For the top mass we have chosen three values which bracket the allowed 
range of variation, namely Mt = 150 GeV, Mt = 175 GeV and Mt = 200 GeV. In table 1 
we report the value of the total cross section, separating the contribution of the gg and 
qq channels. In the first case we further differentiate between the double resonant 
contribution, corresponding to the production of two nearly on-shell top quarks, and 
the background which is dominated by events in which only one top is produced close 
to its mass shell. In the qq channel we distinguish the standard mechanism in which 
the initial pair of light quarks annihilates to a gluon from the O(a^) contribution in 
which the qq pair annihilates to a photon or a Z boson. In both gg-initiated processes 
all non double resonant contributions are negligible. The rightmost column gives the 
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sum of all partial results. All contributions not mentioned in table 1 are negligible, well 
below the percent level. 

Table 1 shows, as well known, that qq annihilation is the dominant source of 
bbW~^W~ events at the Tevatron, while the gg fusion channel contributes about 20% for 
Mt — 150 GeV and only about 7% for Mt — 200 GeV, due to the softer gluon spectrum 
at these energies. At the LHC this pattern is reversed, the gluon fusion mechanism 
provides the bulk of all events while the qq channel contributes between 10 and 20% 
depending on the top mass. We notice that the O(a^) qq annihilation cross section is 
about 2.5% of the corresponding C(q;^q;^) results, independently of the top mass and of 
the coUider energy. The gluon fusion background gives an important contribution to the 
cross section for gg — > bbW'^W~ at the Tevatron, between 6% for Mt — 150 GeV and 
18% for Mt = 200 GcV. However it only contributes about 1% to the total cross section 
for tt production. At the LHC the gg background is smaller compared to the gg signal 
than at the Tevatron, but its contribution to the total ti cross section is larger. This 
can be appreciated in table 2 where we present the contribution in percent to the total 
cross section at the LHC of the three main subprocesses, the double resonant channel 
gg ^ tt ^ bbW~^W~, the single resonant channel gg tbW~ ,ibW~^ — > bbW''^W~ and 
the double resonant qq contribution qq ^ ti bbW^W^ , for the two energy values and 
the three values of the top mass which we have studied. The gluon fusion background 
contributes between 4% and 6% of the total gg cross section, the largest contribution 
being for the heaviest top mass. Comparing the second and third column in table 2 we 
see that the gg background cross section is about one half of the total qq annihilation 
cross section. 

In fig. 3,4 and 5 we present the transverse momentum distribution of the 6-quark 
for Mt = 175, 150 and 200 GeV, respectively at = 10 and 14 TeV. In fig. 6 
one can find the transverse momentum distribution of the for Mt — 175 GeV 
at both LHC energies. The two double resonant contributions and the gluon-gluon 
fusion background are shown separately and the full line gives the sum of the three 
distributions. Within the statistical errors of our Montecarlo the distributions of the 
two double resonant channels differ only in overall normalization and have the same 
shape. The gg background spectrum is peaked at small transverse momenta and 
shows a small hump at about 50 GeV. This can be easily understood noticing that this 
channel is dominated by diagrams, for instance diagram 6 and 7 in fig. 1, in which the 
two initial state gluons split into bb pairs, producing fe's with small Pt- A b from one 
pair and a b from the other then further interact annihilating to a W~^W~ pair through 
the exchange of a ^-channel top. The large mass of the top prevents any substantial 
enhancement for small angle scattering and indeed the background pj^ spectrum does 
not show any preference for small transverse momenta. It is however softer than the 
corresponding double resonant distribution, as can be seen if the different contributions 
are scaled to a common height, say, of the maximum. The hump at large p^ is due to 
the presence of a single top in most of the background events. 

A more refined analysis of the influence of backgrounds on tt production would 
require the introduction of a realistic set of experimental cuts. However, as already 
mentioned, the vast majority of background events contain one nearly on shell t-quark 
and they would survive most selection procedures optimized for the observation of top. 
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leaving the signal to background ratio practically unchanged. A possible catch in the 
preceding argument has to do with the fact that the 6-quark which is not produced in 
the decay of a top, has a soft pt spectrum and some of the tracks from its own decay 
might miss the vertex detector, which is the most effective tool for 6-tagging, whose 
geometrical coverage is hmited. This might reduce the 6-tagging efficiency for this 
channel in comparison with the double resonant one where both 6's are equally hard. 
It is however extremely difficult to estimate this effect without a detailed simulation of 
b hadronization and of the full detector. 

Conclusions 

We have produced the complete matrix elements for gg,qq bbW~^W~. With them 
we have analyzed tt production at the Tevatron and the LflC. With our approach 
we avoid separating the production process from the subsequent decay. Finite width 
effects, irreducible backgrounds and correlations between the two t-quark decays are 
included in our treatment. 

We found that contributions from irreducible backgrounds provide about 1% of the 
total ti cross section at the Tevatron and become more relevant at the LHC, where 
they amount to about 5% of the total production rate. The reason for this difference 
is that single resonant channels are important for the gg cross section and not for the 
qq processes. The the C(q;^) qq annihilation cross section is about 2.5% of the total at 
the Tevatron and is neghgible at the LHC. 
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Table 1 Tree level total cross sections in pb for pp — > bbW'^W at ^/s = 1.8 TeV, 

pp bbW+W~ at ^ = 10 TcV and pp bbW+W' at = 14 TcV with 
^top = 150, 175,200 GeV. The first and second column refer to the contribution 
of double resonant diagrams and to the contribution of single resonant diagrams 
and of the interference between single and double resonant diagrams in gluon- 
gluon fusion respectively. The third and fourth column report the contribution 
from double resonant production via qq fusion into a gluon and into 7 and Z° 
bosons respectively. The last column gives the sum of the four contribution. 

Table 2 Tree level total cross section in pb and contribution in percentage of the most 
relevant channels to the total cross section for pp — > bbW^W~ at ^/s = 10 TeV 
and at ^/s = 14 TeV with m^p = 150, 175, 200 GeV. The first contribution arises 
from gluon fusion double resonant diagrams. The second one is the contribution 
from single resonant diagrams and from the interference between single and double 
resonant diagrams in the gluon fusion channel, while the third one derives from 
quark fusion double resonant diagrams. All other contributions are below 1% at 
tree level. 
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Figure Captions 

Fig. 1 Sample Fcynman diagrams contributing in the lowest order to gg — > bbW'^W~ . 
External wavy lines represent as appropriate. Internal wavy lines represent 
and 7. Gluons permutation are not shown. In diagrams (3)- (5) only t-rcsonant 
diagrams are drawn: t-rcsonant diagrams are obtained replacing b and ^ 
W~ . All other non-resonant diagrams can be obtained from (7) -(8) changing the 
insertions of the internal vector boson hues {Zq and 7). SM Higgs contributions 
are not shown because their contribution is far below the accuracy reached in the 
present work. They can be obtained replacing internal neutral vector boson lines 
with scalar Higgs lines. 

Fig. 2 Sample Feynman diagrams contributing in the lowest order to qq — > 661^ "'"1^^. 
External wavy lines represent as appropriate. Internal wavy lines represent 
Z° and 7. Gluons permutation are not shown, f-resonant diagrams are obtained 
replacing 6 6 and ^ W~ in diagram (2). All other non- resonant diagrams 
can be obtained from (3) -(6) changing the insertions of the internal vector boson 
lines (Zo and 7). SM Higgs contributions are not shown because their contribution 
is far below the accuracy reached in the present work. They can be obtained 
replacing internal neutral vector boson lines with scalar Higgs lines. 

Fig. 3 The distribution of transverse momentum of the b quark for pp — > bbW'^W~ 
with mtop = 175 GeV at = 10 TeV (Fig. 3a) and ^ = 14 TeV (Fig. 3b). 
The dotted line represents the contribution from double resonant diagrams in 
the quark fusion channel, while the dash-dotted line is the contribution from 
double resonant diagrams in gluon fusion. The dashed line is the "irreducible 
background" from single resonant diagrams and from the interference between 
single and double resonant diagrams in the gluon fusion channel. The continuous 
line is the sum of the three contributions. 

Fig. 4 The distribution of transverse momentum of the b quark for pp — > bbW'^W~ 
with mtap = 150 GeV at = 10 TeV (Fig. 4a) and = 14 TeV (Fig. 4b). 
Symbols as in Fig. 3. 

Fig. 5 The distribution of transverse momentum of the b quark for pp bbW^W~ 
with m^p = 200 GeV at = 10 TeV (Fig. 5a) and = 14 TeV (Fig. 5b). 
Symbols as in Fig. 3. 

Fig. 6 The distribution of transverse momentum of the for pp bbW^W^ with 
mtop = 175 GeV at ^ = 10 TeV (Fig. 6a) and y/s ^ U TeV (Fig. 6b). Symbols 
as in Fig. 3. 
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Fig.2 
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